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Abstract 

In this paper it is investigated how to find a matrix representation of operators on 
a Hilbert space Ti with Bessel sequences, frames and Riesz bases. In many applica- 
tions these sequences are often preferable to orthonormal bases (ONBs). Therefore 
it is useful to extend the known method of matrix representation by using these 
sequences instead of ONBs for these application areas. We will give basic defini- 
tions of the functions connecting infinite matrices defining bounded operators on I 2 
and operators on TC. We will show some structural results and give some examples. 
Furthermore in the case of Riesz bases we prove that those functions are isomor- 
phisms. Finally we are going to apply this idea to the connection of Hilbert-Schmidt 
operators and Frobenius matrices. 



1 Introduction 



The relevance of signal processing in today's life is clearly evident. Without 
exaggeration it can be said, that any advance in signal processing sciences 
directly leads to an application in technology and information processing. 
Without signal processing methods several modern technologies would not be 
possible, like mobile phone, UMTS, xDSL or digital television. 

The mathematical background for today's signal processing applications are 
Gabor [15] , wavelet p2] and sampling theory [6]. A signal is sampled and then 
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analyzed using a Gabor respectively wavelet system. Many applications use a 
modification on the coefficients obtained from the analysis operation [20f2~T] . 
For them not only an analysis but also a synthesis operation is needed. If the 
coefficients are not changed, the result of the synthesis should be the original 
signal, therefore so-called perfect reconstruction is needed. One way to achieve 
that is to analyze the signal using orthonormal bases (ONBs). In this case 
the analysis of a function is simply the correlation of the signal / with each 
basis element e&, / i— > (c*.) := ((/, e&)) . The synthesis that gives perfect 
reconstruction is simply the (possibly infinite) linear combination of the basis 
elements using the coefficients c = (c^), ch>^ c^k- 



From practical experience it soon has become apparent that the concept of 
an orthonormal basis is not always useful. Sometimes it is more important for 
a decomposing set to have other special properties rather than guaranteeing 
unique coefficients. For example it is impossible to have good time-frequency 
localization for Gabor ONBs or a wavelet ONB with a mother wavelet which 
has exponentially decay and is infinitely often differentiable with bounded 
derivatives [10j . Furthermore suitable ONBs are often difficult to construct 
in a numerical efficient way. This led to the concept of frames, which was 
introduced by Duffin and Schaefer [T5] . It was made popular by Daubechies 
[T2] . and today it is one of the most important foundations of Gabor [T5] . 
wavelet [2] and sampling theory [TJ. In signal processing applications frames 
have received more and more attention [71128] . 

Models in physics [2] and other application areas, for example in sound vi- 
bration analysis [5], are mostly continuous models. A lot of problems there 
can be formulated as operator theory problems, for example in differential or 
integral equations. To be able to work numerically the operators have to be 
discretized. One way to do this is to find (possibly infinite) matrices describ- 
ing these operators using ONBs. In this paper we will investigate a way to 
describe an operator as a matrix using frames. This kind of 'sampling of op- 
erators' (compare to [23]) is especially important for application areas, where 
frames are heavily used, so that the link between model and discretization is 
kept. For implementations operator equations can be transformed in a finite, 
discrete problem with the finite section method [T7J in the same way as in the 
ONB case. 

The standard matrix description [TTJ of operators O using an ONB (e^) is by 
constructing an matrix M with the entries Mj^ = (Oek,ej). In [9] a concept 
was presented, where an operator R is described by the matrix ((R<pj, (pi)) . 



with (<f>i) being a frame and ((pi) its canonical dual. Recently such a kind 
of representation is used for the description of operators in [19] using Gabor 
frames and [25] using linear independent Gabor systems. In this paper we are 
going to develop this idea in full generality for Bessel sequences, frames and 
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Riesz sequences and also look at the dual function which assigns an operator 
to a matrix. 

This paper is organized as follows: In Section [2] we collect results and nota- 
tion we need. Section 13.11 gives the basic definitions and properties for Bessel 
sequences and frames. Matrix representation with Riesz bases is covered in 
Section 13.21 In Section 13.31 the connection of Frobenius matrices and Hilbert- 
Schmidt is investigated. Section @] finishes the paper with perspectives. 



2 Notation and Preliminaries 

2.1 Hilbert spaces and Operators 

We will give only a short review, for details refer to [11] . We will denote infinite 
dimensional Hilbert spaces by Ti and their inner product with < ., . >, which 
is linear in the first coordinate. Let i3(7^i,7^2) denote the set of all linear 
and bounded operators from TC\ to Ti.2- With the operator norm, \\A\\ = 

sup | ||j4(x) \\ H2 j, this set is a Banach space. We will denote the composition 

of two operators A : Hi — > H2 and B : Ti 2 — > TC3 by B o A : TCi — > H3 and the 
adjoint of an operator A by A*, so that (Ax, y) = (x, A*y) for all x, y G Ji. 

Furthermore we will denote the range of an operator A by ran(0) and its 
kernel by ker(A). An example for a Hilbert space is the sequence space I 2 
consisting of all square-summable sequences in C with the inner product 
(c,d) = Y^ c k ■ dk- We will use the canonical basis A = (5k) for sequence 

k 

!1 k = n 
otherwise 

Remember that a linear function between Banach algebras tp : Q3i — > 2^2 is 
called a Banach algebra homomorphism, if it is also multiplicative, i.e. for all 
x, y G Q3i we have tp(x ■ y) = tp(x) ■ tp(y). It is called a monomorphism, if it is 
also injective. 

Definition 2.1 Let X,Y,Z be sets, f : X — > Z , g : Y — ► Z be arbitrary 
functions. The Kronecker product ® : X x Y — > Z is defined by 

(f®o 9) (x,y) = f(x) ■ g(y). 
Let f G TCi, g G TC2 then define the inner tensor product as an operator from 
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Ti.2 to Hi by 

(f® t g) (h) = (h,g)f forheH 2 . 



We will often write / <S> g instead of / ® Q g or / <g>j g , if there is no chance of 
misint erpret at ion . 



2.1.1 Hilbert Schmidt Operators 

A bounded operator T e B(H\,H2) is called a Hilbert- Schmidt (HS) operator 
if there exists an ONB (e n ) C 7Yi such that 



ITII 



\ £ W Te n\\n 2 < 00 

\ n=l 



Let HS(H\,H2) denote the space of Hilbert Schmidt operators from H\ to 
H 2 . 

This definition is independent of the choice of the ONB. The class of Hilbert- 
Schmidt operators is a Hilbert space of the compact operators with the fol- 
lowing properties: 

• ll^llop — Ik II-hs 

• \\T\\ HS = \\T*\\ HS , and T e HS ^ T* e HS. 

• If T E HS and A e B, then TA and AT e HS. \\AT\\ HS < \\A\\ 0p \\T\\ ns 
and \\TA\\ HS < ||A|| 0p ||r|| ws . 

For more details on this class of compact operators refer to [23] or [29] . 



2.2 Frames 



For more details and proofs for this section refer e.g. to [8|10f 121118] . 

A sequence \& = (tpk\k £ K) is called a frame for the Hilbert space H, if 
constants A, B > exist, such that 

A-\\f\\ 2 H <T,\(fM\ 2 < B -Wf\\H V/eH (1) 



Here A is called a lower, B an upper frame bound. If the bounds can be chosen 
such that A = B the frame is called tight. 

A sequence ^ = (Jj^) is called a Bessel sequence with Bessel bound if it 
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fulfills the right inequality above: 



J2\(f^ k )\ 2 <B-\\f\\ 2 H VfeH (2) 

k 



The index set will be omitted in the following, if no distinction is necessary. 

For a Bessel sequence, \& = (ipk), let C-q, : 7i — > 1 2 {K) be the analysis operator 

C*(/) — ((/)^fc))fc- Let : Z 2 (i^) — > 7ibethe synthesis operator Dq, ((c k )) = 

Ecfc • V'fc- Let : — > be the (associated) frame operator Sq,(f) = 
k 

E (/) V'fc) " V'fe- To simplify notation we will just write S for Sq?, C for and 

k 

D for D^, if it is not necessary to distinguish different frames. We will use the 
notation Sqr : q, = Dq, o C$. C and D are adjoint to each other, D — C* with 
1 1 -^11 o P = W^Wop — The series Sc fc - ip k converges unconditionally for all 

(c k ) e I 2 . 

For a frame ^ = (^) with bounds A, C is a bounded, injective operator 

with closed range and S = C*C = DD* is a positive invertible operator 

satisfying AI n < S < BI n and B~ l I n < S^ 1 < Even more we can 

find an expansion for every member of H: The sequence \I> = (ipk) = (S^ipk) 

is a frame with frame bounds B^ 1 , A^ 1 > 0, the so called canonical dual frame. 

Every f E H has the expansions / = E (f,i>k)i>k and / = E {f,*Pk)4>k 

k&K x 7 feeK 
where both sums converge unconditionally in 7Y. 

Remember that a sequence (e^) is called a (Schauder) basis for 7i, if for all 

f E H there are unique coefficients (c k ) such that / = J2 c k<Pk- Also two 

fe 

sequences (ip k ), (4>k) are called biorthogonal if (ip k ,<pj) = Skj f° r all 

A complete sequence (^) in 7i is called a Riesz basis if there exist constants 
A, £> > such that the inequalities 



^IM|2< 



c ^ 

keK 



<B\\c\\ 2 2 



H 



hold for all finite sequences (c k ). 



For a frame (ip k ) the following conditions are equivalent: (i) (ip k ) is a Riesz 
basis for H. (ii) The coefficients (c k ) G / 2 for the series expansion with (ip k ) are 
unique. So the synthesis operator D is injective. (m) The analysis operator C 
is surjective. (iv) (ipk) and (^fc) are biorthogonal. 

Let ^ = (ipk) and $ = be two sequences in Ti. The Gram matrix Gq,^ for 
these sequences is given by (G^,$) jm = (4> m ,ipj), j,m G if. We denote 
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by G^. We can look at the operator induced by the Gram matrix, defined for 
eel 2 formally as (G*,$c)j = J2 c k {4>k,ipj) ■ Clearly for two Bessel sequences 

k 

it is well defined as linear bounded operator, because 

(Gvpcy = Y c k ^j) = ( Y ipj ) = {{C 9 o D$) c) j 

k \ k I 

and therefore ||G*,$|| 0p < ||C*||o p ||-^*Hop — A frame is a Riesz sequence 
if and only if the Gram matrix defines a bounded and invertible operator on 
I 2 . 



3 Representing Operators with Frames 

Let (ipk) be a frame in TCi- An existing operator U E B(TCi,H.2) is uniquely 
determined by its images of the frame elements. For f = Ck*Pk 

k 

u(f) = uC£cM = Y c kU^k). 

k k 

On the other hand, contrary to the case for ONBs, we cannot just choose a 
Bessel sequence (rjk) and define an operator just by choosing V(V'fc) := Vk and 
setting V(Y, Ckipk) — Y^ c kVk- This is in general not well-defined. Only if 

k k 

Y C ki>k = Y d ^k Y C kVk = Y d kVk 

k k k k 

this definition is non-ambiguous, i.e. if ker (D^ k ) C ker (D Vk ). This condition 
is certainly fulfilled, if D^ k is injective, i.e. for Riesz bases. 

This problem can be avoided by using the following definition 

V(f):=Y(f^k)vk. (3) 

k 

As (rjk) forms a Bessel sequence, the right hand side of Eq. d3J) is well-defined. 
It is clearly linear, and it is bounded. The Bessel condition, Eq. is necessary 
in the case of ONBs to get a bounded operator, too [UJ. But contrary to the 
ONB case, here, in general, V(ipk) ^ Vk- 

Instead of changing the sequence with which the coefficients are resynthezised, 
an operator can also be described by changing the coefficients, as presented 
in the following sections. 
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3.1 Matrix Representation 



For orthonormal sequence it is well known, that operators can be uniquely 
described by a matrix representation [17]. The same can be constructed with 
frames and their duals. Recall the definition of the operator defined by a 
(possibly infinite) matrix: (Mc) , = SM^c^. We will start with the more 

k 

general case of Bessel sequences. Note that we will use the notation H-H^^.^ 
for the operator norm in B(Ti.i,Ti.2) to be able to distinguish between different 
operator norms. 

Theorem 3.1 Let \& = (ip^) be a Bessel sequence in TCi with bound B, $ = 
in Ti.2 with B' . 

(1) Let O : 7i\ — > Ti.2 be a bounded, linear operator. Then the infinite matrix 

defines a bounded operator from I 2 to I 2 with ||A^||;2_^.,2 < V B ■ B' ■ 

M (*,v) (O) = o O o Aj> 

This means the function JA^' 9 ' : Bi^Hx.H.^) — > B(l 2 ,l 2 ) is a well-defined 
bounded operator. 

(2) On the other hand let M be an infinite matrix defining a bounded operator 
from I 2 to I 2 , (Mc). = X) MfcCfc- Then the operator 0(*>*) defined by 

k 

( M )) & = £ fe M kJ (h, V>J> \ <Pk, for heHr 

is a bounded operator from Tlx to 7^2 with 

0<*'*>(M) < VbTW' \\M\\ P ^ P . 

Hi — >rL2 

£>(*>*) (M) = D$oMoC* = J2Y, M k,j ■ <Pk ®* % 

k j 

This means the function 0(*>*) : B(l 2 ,l 2 ) — > B(TCi,Ti.2) is a well-defined 
bounded operator. 

Proof: Let M = and O = 0^>*\ Let O G B(Hi,H 2 ), then 

(M (O) C ) =Y,(M (0)) jtk c k = £ <OVfc, 0i> = 
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O C k^k 
\ k 



(4) 



Il-^ c ll2 = E l(0^c, 0,,->| 2 < B' • ||OI>*c|£ < 5' ■ ||0||^ B 

j 

Equation H] also shows us, that as operator we have 

M (*,9) (0) = C$oOoDy. 



On the other hand let M be an infinite matrix, then 

O (M) = D$ o M o 

<v^-||m|| m Vb 

Hi ^ W a 
M *" 



□ 



Fig. 1. The operator induced by a matrix M and the matrix induced by an operator 
O 



Definition 3.1 For an operator O and a matrix M as in Theorem lff.il we 
call My 5, \0) the matrix induced by the operator O with respect to the Bessel 
sequences \I> = (/0fc) and $ = and 0^>®>{M) the operator induced by the 
matrix M with respect to the Bessel sequences \& and $. (See Figure^) 

If we do not want to stress the dependency on the frames and there is no 
change of confusion, the notation M.(0) and 0(M) will be used. 

For frames we can prove more properties: 

Proposition 3.2 Let ^ = (i/j^) be a frame in TCi with bounds A, B, $ = (0^) 
in H.2 with A',B' . Then 

(1 ) o M@>*A =Id= o Af . 

ylnd therefore for all O E B(TCi,TC2}-' 



o = J2 (oii>j, <t>k) 4>k ®i ipj 

k,j 

(2) Ai^'^ is infective and 0(*>*) is surjective. 

(3) Let H 1 = H 2 , then 0(*'*)(ld {2 ) = 
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(4) Let E = be any frame in H3, and O : H3 — > Hi and P : Hi — > H3. 
Then 

M (*M {OoP)= (M ( *' s) (O) ■ (P)) 
Proof: 1.) For / G Hi we get (O&M o M^) (O) (f) = E ( E (C>4 ~4>k) (f, ^) 

& \ 3 

= E ((e (/. ^> o^, ^ j fc = E (o/, 0*) 0^ = o/. 

For the other equality the roles of the frame and the dual just have to be 
switched. 

2. ) From OM. = Id we know that M. is injective and O is surjective. 

3. ) 0(ld)f = E f E (/, i)) ^fc = E (/, ^ fc ) V* = /• 

4. ) (O o p) m = (O o P^ ?! p > = (P^ ?! 
On the other hand 

(a*<*- s > (o) ■ (p)) = E-^ ($,H) (o) p . fc • jw^*) (p) m = 
= 2 (o/ifc, P ) (py>„ ik) = E (p^ ik) = 

k k 




□ 

As a direct consequence we get the following corollary: 

Corollary 3.3 For the frame $ = the function JA^'^ is a Banach- 

algebra monomorphism between the algebra of bounded operators (B(Hi,Hi), °) 
and the infinite matrices of {B(l 2 , I 2 ), ■). 

The other function O is in general not so "well-behaved". It is, if the dual 
frames are biorthogonal. In this case these functions are isomorphisms, refer 
to Section 13.21 

Lemma 3.4 Let O : Hi — > H2 be a linear and bounded operator, let = (ip k ) 
and $ = (0^) be frames in Hi resp. Hi- Then AA^'^\0) maps ran (CV) into 
ran (C$) with 

{{fM) k ^{{OfAk)) k - 



9 



If O is surjective, then Ai^'^(O) maps ran{C\&) onto ran(C§). If O is in- 
jective, M. (O) is also injective. 

Proof: Let c G ran{Cq), then there exists / G TCi such that Ck = (f,ipk)- 

(m^\o)(c)). = J2(o$ k ,<f>i) (fM = (e(/>^)°^) = (of,<j>i) 

k ' \ k I 

So ((f^k)) k ^((Of,cp k )) k . 

If O is surjective, then for every / there exists a g such that Og = f, and 
therefore (g,ip k ) ^ (f,<Pk)- 

If is injective, then let's suppose that (Of,<pk) = (Og,4>k)- Because {<pk) is 
a frame => Of = Og => f = g =>• (f,ip k ) = (g, ip k )- D 



Particularly for O = Id the Gram matrix G^ ^ = ((ip k , </>i)j maps ran (C#) 

bijectively on ran (C$). So we get a way to a way to "switch" between frames 
by mapping from one analysis range into the other [3] . 

Let us give some examples: 

Example 3.1 : 

Let \I/ = {ipk) and $ = {<p k ) be frames in and A = {5 k ) the canonical basis 
of I 2 . Then 

(1) S^-.n^H^ M = 

(2) S^ 1 : w - n =► .M^W) = G*. 

(3) C$ : ft -> Z 2 

^^(C*)^ = (CftV'i, 4) = E frMi) (4, S t ) = fa, 4> k ) = (G*,*) 



(4) id-.n^n^ MW(i<t) = g 9>9 . 

(5) Id:l 2 ^l 2 ^ (Id) = D$ o C § = ^. 



5. i.i Motivation: Solving Operator Equalities 
Given an operator equality 

0-f = g (5) 

it is natural to discretize it to find a solution. Let $ = be a frame. Let us 
suppose that for a given g with coefficients d = {d k ) = {{g, 4>k)) an d a matrix 
representation M of O there is an algorithm to find the least square solution 
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of 



M-c = d (6) 

for example using the pseudoinverse [10]. Still, if using frames, we can not 
expect to find a true solution for Eq. [5] just by applying D^, on c as in general 
c is not in ran(C$) even if d is. But rephrasing Eq. [5] we see the following: 

Of = g<^> y E(f> °^ = 9 E </> = 0*> 

fc fc 

<=► (O) ■ C./ = C 9 g. 

It can be easily seen that this is equivalent to projecting c on ran(C), solving 
MC<t>Dq,c = d, which is a common idea found in many algorithms, for example 
for a recent one see [26J. 

This gives us an algorithm for finding an approximative solution to the inverse 
operator problem Of = g. 

(1) Set M = M^(0). 

(2) Find a good finite dimensional approximation Mn of M by using the 
finite section method [T7] and 

(3) then apply an algorithm like e.g. the QR factorization [27] to find a 
solution for Eq. El 

(4) and synthezise with the dual frame $. 

Remark: It has been shown in [22], that the finite section is very useful in 
the case of frame theory. It would be very interesting to investigate the idea 
presented above further in this context. 



3.2 Matrix representation using Riesz Bases 

The coefficients using a Riesz basis are unique, so Theorem 13 . 1 1 can be extended 
to: 

Theorem 3.5 Let $ = be a Riesz basis for Ttx, \P = (ipk) one for Ti.2- 
The functions Ai^'®' and between B(TCi, and the infinite matrices 

in B(l 2 ,l 2 ) are bijective. JA^'^ and are inverse to each other. For 

7i\ = 7^2 the identity is mapped on the_ identity by Ai^'^ and If 
furthermore \I> = $ then Ai^'^ and C?(*>*) are Banach algebra isomorphisms, 
respecting the identities idp and idyi- 
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Proof: We know from Proposition 13.21 that O o M. = Id. Let's look at 

((M o O) {M)) M = M [EE M Ki (; 1>i) n 



p,q 



E E Mfc j ^t.^ =EE M k,i (&> P ) = M P,g 

\ k j I k j 



So these functions are inverse to each other and therefore bijective. 



We know that A4^'^ is a Banach algebra homomorphism and so is its inverse. 
□ 



3.3 Matrix Representation ofTiS Operators 

We now have the adequate tools to state that HS operators correspond exactly 
to the Frobenius matrices, as expected. 

Definition 3.2 Let A be an m by n matrix, then 

\\A\ 



'fro 



n—1 m—1 



y E E 

\ i=0 j=0 



is t/ie Frobenius norm. Let us denote the set of all matrices with finite Frobe- 
nius norm by l^ 2 ' 2 ', the set of Frobenius matrices. 

Proposition 3.6 Let ^ = (ipk) be a Bessel sequence in Hi with bound B, 
$ = (0 fc ) i n n 2 with B'. Let M be a matrix in l (2 > 2) . Then 0^{M) e 
HS (7ix, H-i), the Hilbert Schmidt class of operators from Hi to H2, with 
\\0(M)\\ ns < VBBl\\M\\ fro . 

Let O e HS, then M&> 9 \0) G lS 2 ^ with \\M(0)\\ fro < VBB 1 \\0\\ ns . 

Proof: 1.) Naturally the matrices in l^ 2,2 ^ correspond to Hilbert-Schmidt 
operators on I 2 as \\M\\ 2 HS(l2] = E HAf^l* = EE = ||M|L . As the 

' i i P 

Hilbert-Schmidt class of operators is an ideal, we know that 
\\0(M)\\ ns{HliH2) = \\D* o M o C*\\ ns{Hi>H2) < 



12 



< \\D*\\op(P,n>) ■ W M Whs(i>) ■ IICWllppCKi^) = ^BB'\\M\\ fro 
2.) Let (e p ) be any ONB of H 2 . 

\\M(o)f fro = EE iW^>l 2 <E 5 '- \\oM 2 ni = E^'EIWWI 2 = 

ik i ip 

= E £'EK^ 0%>| 2 < E HO*e p ||^ 2 < B'B ||0*||^ = Bff \\0\\ 2 ns 
v i v 

□ 



3.3.1 Matrices and the Kernel Theorems 



For L 2 {W l ) the TCS operators are exactly those integral operators with kernels 
in L 2 (]R M ) [HEM] . This means that there exists a k G L 2 (IR 2d ) such an 
operator can be described as 

(Of) 0) = / K (x,y)f(y)dy 

Or in weak formulation 

(Of, 9) = J J ko(x, y)f\y)g{x)dydx = (k , f ® g) . (7) 
From [37X1 we know that 

o = E ( $?> 0*) 0* ®< 

and so 

Corollary 3.7 Let O eUS (L 2 (R d )) . Let * = (^) and $ = (<f> k ) be frames 
in L 2 (M. d ^j . Then the kernel of O is given as: 

K o = E fa) ' fa ®o vpj = E " fa ®o iij 

j,k j,k 

Proof: 

K(0) = K E 0fc) ^ ®i 

= E 0*) K ®< ^j) = E (O^jifa) <Pk Oo 

as (/ ®i g) (h) = (h, g) f = f h(x)g(x)dxf(y) and so «(/ ® t g) = f <S> 

□ 



13 



4 Perspectives 



Apart from TiS operators on L 2 (M. d ), there is a large variety of function spaces, 
where operators are exactly integral operator using the weak formulation in 
Eq. ([7]), for example for operators on the Schwartz space O : S — > S' , on 
modulation spaces O : M*(IR d ) — > My v [18], or on Feichtinger's algebra O : 
S' — > So and O : So — > S"o [E] . In order to derive similar results for the case of 
Banach spaces of functions or distributions, Section 13.11 could be generalized 
to these function spaces. 

Using the connection between operators and frames, we can ask, which oper- 
ators are induced by diagonal matrices. Let m be a sequence and Diag(m) the 
matrix that has this sequence as diagonal. Then define 

M mi $,f := C) ( *'* ) (Diafl(m)) = X) m * ' 0* ® V'fe 

This means we have arrived quite naturally at the definition of frame multiplier 
as introduced in [I]. This connection should be investigated and exploited 
further. 

Especially in areas where orthonormal bases are not very useful and frames 
are already used heavily for analysis-synthesis systems , the idea of matrix rep- 
resentation with frames should be used in applications. This concept should 
be applied using the finite section method [17] as mentioned in Section 13.1.11 
Furthermore, to conserve the semantic connection between model and dis- 
cretization even more, it might be interesting to use frames in a projection 
method like the finite section method for approximation of infinite matrices 
by finite ones. Instead of using the ONB for the projection use the frames, 
which are already connected to the problem. In this sense the ideas in [22] 
could be extended to the notion of self-localized frames as presented in [16J. 



5 ACKNOWLEDGEMENT 



The author would like to thank Hans G. Feichtinger, Bruno Torresani and 
Jean-Pierre Antoine for many helpful comments and suggestions as well as 
Wolfgang Kreuzer for proofreading. He would like to thank the hospitality 
of the LATP, CMI, Marseille, France and FYMA, UCL, Louvain-la-Neuve, 
Belgium, where part of this work was prepared, supported by the HASSIP- 
network. This work was partly supported by the European Union's Human 
Potential Program, under contract HPRN-CT-2002-00285 (HASSIP). 



14 



References 



[I] A. Aldroubi and K. Grochenig. Non-uniform sampling and reconstruction in 
shift-invariant spaces. SIAM Review, 43:585-620, 2001. 

[2] S. T. Ali, J. -P. Antoine, and J. -P. Gazeau. Coherent States, Wavelets and Their 
Generalization. Graduate Texts in Contemporary Physics. Springer New York, 
2000. 

[3] P. Balazs. Regular and Irregular Gabor Multipliers with Application to 
Psychoacoustic Masking. Phd thesis, University of Vienna, June 2005. 

[4] P. Balazs. Basic definition and properties of Bessel multipliers. Journal of 
Mathematical Analysis and Applications, 325(l):571-585, January 2007. 

[5] P. Balazs, W. Kreuzer, and H. Waubke. A stochastic 2d-model for calculating 
vibrations in liquids and soils, accepted for Journal of Computational Acoustics, 
2006. 

[6] J. Benedetto and P. Ferreira. Modern sampling theory. Mathematics and 
applications. Birkhauser, 2001. 

[7] H. Bolcskei, F. Hlawatsch, and H. G. Feichtinger. Frame-theoretic analysis of 
oversampled filter banks. IEEE Trans. Signal Processing, 46(12) :3256-3268, 
1998. 

[8] P. G. Casazza. The art of frame theory. Taiwanese J. Math., 4(2):129-202, 
2000. 

[9] O. Christensen. Frames and pseudo-inverses. J. Math. Anal. Appl, 195(2):401- 
414, 1995. 

[10] O. Christensen. An Introduction To Frames And Riesz Bases. Birkhauser, 2003. 

[II] J. B. Conway. A Course in Functional Analysis. Graduate Texts in 
Mathematics. Springer New York, 2. edition, 1990. 

[12] I. Daubechies. Ten Lectures On Wavelets. CBMS-NSF Regional Conference 
Series in Applied Mathematics. SIAM Philadelphia, 1992. 

[13] R. J. Duffin and A. C. Schaeffer. A class of nonharmonic Fourier series. Trans. 
Amer. Math. Soc, 72:341-366, 1952. 

[14] H. G. Feichtinger and W. Kozek. Quantization of TF lattice-invariant operators 
on elementary LCA groups, chapter 7, pages 233-266. In [15], 1998. 

[15] H. G. Feichtinger and T. Strohmer. Gabor Analysis and Algorithms - Theory 
and Applications. Birkhauser Boston, 1998. 

[16] M. Fornasier and K. Grochenig. Intrinsic localization of frames. Constructive 
Approximation, 22(3):395-415, 2005. 



15 



[17] I. Gohberg, S. Goldberg, and M. Kaashoek. Basic Classes of Linear Operators. 
Birkhauser, 2003. 

[18] K. Grochenig. Foundations of Time- Frequency Analysis. Birkhauser Boston, 
2001. 

[19] K. Grochenig. Time-frequency analysis of Sjostrand's class. Rev. Mat. Iberoam., 
22: (to appear), 2006. 

[20] R. Kronland-Martinet and A. Grossmann. Application of time-frequency 
and time-scale methods (wavelet transform) to the analysis, synthesis and 
transformation of natural sounds. In G. De Poli, A. Piccialli, and C. Roads, 
editors, Representations of musical signals, pages 45-85. MIT Press, Cambridge, 
1991. 

[21] G. Matz and F. Hlawatsch. Linear Time- Frequency Filters: On-line 
Algorithms and Applications, chapter 6 in 'Application in Time-Frequency 
Signal Processing', pages 205-271. eds. A. Papandreou-Suppappola, Boca 
Raton (FL): CRC Press, 2002. 

[22] O.Christensen and T.Strohmer. The finite section method and problems in 
frame theory. Journal of Approximation Theory, 133(2):221-237, 2005. 

[23] G. Pfander and D. Walnut. Operator identification and Feichtinger's algebra, 
accepted for Sampling Theory in Signal and Image Processing, 2006. 

[24] R. Schatten. Norm Ideals of Completely Continious Operators. Springer Berlin, 
1960. 

[25] T. Strohmer. Pseudodifferential operators and Banach algebras in mobile 
communications. Appl.Comp. Harm. Anal., in press. 

[26] G. Tcschke. Multi-frame representations in linear inverse problems with mixed 
multi-constraints. DFG-SPP-1114 preprint 90, (submitted), 205. 

[27] L. N. Trefethen and D. Bau III. Numerical Linear Algebra. SIAM Philadelphia, 
1997. 

[28] M. Vetterli and J. Kovacevic. Wavelets and Subband Coding. Prentice Hall, 
1995. 

[29] D. Werner. Funktionalanalysis. Springer Berlin, 1995. 



16 



